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Abstract 

We consider Euclidean Conformal Field Theories perturbed by quenched disorder, 
namely by random fluctuations in their couplings. Such theories are relevant for second- 
order phase transitions in the presence of impurities or other forms of disorder. Theories 
with quenched disorder often flow to new fixed points of the renormalization group. We 
begin with disorder in free held theories. Imry and Ma showed that disordered free fields 
can only exist for d > 4. For d > 4 we show that disorder leads to new fixed points which 
are not scale-invariant. We then move on to large-A theories (vector models or gauge 
theories in the ‘t Hooft limit). We compute exactly the beta function for the disorder, and 
the correlation functions of the disordered theory. We generalize the results of Imry and 
Ma by showing that such disordered theories exist only when disorder couples to operators 
of dimension A > d/4. Sometimes the disordered fixed points are not scale-invariant, and 
in other cases they have unconventional dependence on the disorder, including non-trivial 
effects due to irrelevant operators. Holography maps disorder in conformal theories to 
stochastic differential equations in a higher dimensional space. We use this dictionary 
to reproduce our field theory results. We also study the leading 1/N corrections, both 
by held theory methods and by holography. These corrections are particularly important 
when disorder scales with the number of degrees of freedom. 
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1. Introduction and Review 

It is of considerable theoretical and experimental interest to consider Qnantnm Field 
Theory (QFT) with random conplings. By “random couplings” or “quenched disorder” 
we mean that we have some specihed probability distribution in the space of couplings 
constants. In the simplest case one can let the couplings vary independently at every 
point in Euclidean space (though the ensemble should still consist of reasonably smooth 
functions). This probability distribution is independent of the quantum helds, which is 
why this type of disorder is called “quenched.” We imagine hrst solving the theory for all 
values of the couplings, and then averaging over the couplings. 

Experimentally, this situation arises (in some cases) for the conformal held theories 
that describe the critical points of second-order phase transitions, in the presence of disor¬ 
der. Such disorder is often present, due to impurities or to random external helds. Upon 
adding disorder, the phase transition point may be modihed (i.e. the critical exponents 
may change), it can disappear altogether, or it can be unaltered.0 

^ We do not study in this paper disordered theories relevant to quantum phase transitions. 
In these theories there is an extra time direction, but the disorder only depends on the spatial 
directions. We also do not consider thermal phase transitions in quantum systems, which are 


1 





























































Let us first define the problem (see, for instance, Consider a Euclidean QFT 

in d dimensions, in which all the operators appear in the action with couplings that 
are allowed to depend on space, S = Sq-\- f d‘^xg^(x)Oi(x). We can then dehne the 
generating functional W[ 5 r*] for connected diagrams (i.e. minus the free energy) 


eW[9l ^ z[g^] = / [Vg]e 






where [Dg] stands for the path integral measure. In this paper we focus on the case where 
the theory with = 0 is a Conformal Field Theory (CFT). Functional derivatives of IF[ 5 r*] 
evaluated at g'^{x) = 0 give connected correlation functions in the CFT. 

Now, suppose that one of the couplings associated to a scalar operator, h = is 
random, with a Gaussian probability distributioni 


h{x) = 0 , h{x)h{y) = c^d^'^\x — y) , 


( 1 . 2 ) 


where we denote by X the average of X over the disorder. The disorder average of higher 
products of h{xi) is given by a sum over all the two-point functions as in Wick’s theorem. 

We can then dehne the disordered free-energy Wd as an average over the ordinary 
free-energy 

WD[g\g\...] = I , (1.3) 

and the disorder-averaged vacuum-normalized correlation functions are given by 

{Oi{x,)...OM) = e-^ I ‘^"^^'^^\Ol{xi)...Or,{Xr,))H , (1.4) 


where {Oi{xi)...On{xn))h is the usual vacuum-normalized path integral with a source h{x) 


{Oi{xi)...On{Xn))h = 




(1.5) 


The operators Oi may or may not coincide with Oq. 


described by a compact Euclidean time direction. We only consider here the most symmetric case 
of disorder in Euclidean CETs; this is relevant to classical phase transitions. We leave the study 
of the various other cases by our methods to future work. 

^ The case of more than one coupling being random can be treated similarly, as well as cases 
where the disorder couples to operators with non-zero spin. Also, one can consider other proba¬ 
bility distributions in the space of couplings. 
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Of particular interest are the correlation functions that can be obtained by derivatives 
of the disordered free energy, In the context of second order phase transitions these 

correlation functions determine many aspects of the thermodynamics, and thus play an 
important role. Unlike the usual situation in CFTs where all correlation functions can be 
obtained from derivatives of the free energy, this is not the case for disordered theories. 
We will be careful to illuminate these differences when they are important. For example, if 
we take the second derivative of the disordered free energy with respect to some couplings 
( 7 * that couple to Oi we would hnd the disorder-averaged connected correlation function 

_ _ 

.S-(xOVfe) = . (1.6) 


Note that the second term is not the same as {Oi{xi))-{Oj{x 2 ))- Generically the correlation 
functions (|1.6[) behave differently from {Oi(xi)Oj(x 2 ))- In the context of the holographic 
duality of held theories to gravitational theories, connected correlation functions of the 
type ( p.. 6 |) naturally appear as derivatives of the bulk action of the gravitational theory. 

The Harris argument provides a simple diagnostic of when disorder is irrelevant. 
Suppose the scaling dimension of the scalar operator Oq to which we couple disorder is A. 
Because of the quadratic term in h on the right-hand side of (|1.3|) , it is natural to associate 
to h dimension d/2, and to take c to be a dimensionless coupling constant. Then, the term 
f d^xh[x)0{x) in the action is relevant if A < d/2, marginal if A = d/2, and irrelevant 
if A > d/2. Thus, we expect disorder to be innocuous for A > d/2. This is the Harris 
bound. If A < (i/2 we have to analyze the situation case by case. In this situation, even 
arbitrarily small disorder may have profound effects on the low-energy behavior of the 
theory. It is also useful to keep in mind that the argument that for A > d/2 disorder has 
no effect on the phase transition is only valid for small disorder; large disorder is similar 
to an irrelevant operator with a large coefficient, and it may have non-trivial consequences 
(another subtlety in this argument is considered in subsection 2.4). 


One difficulty in analyzing the disordered free energy ( |1.3|) by standard held theory 
methods is that we are averaging the logarithm of the partition function, and not the 
partition function itself, so we are not just adding extra terms to the action. One common 
way to address this is by appealing to the replica trick, where we average the n’th power 
of the partition function, and then perform an analytic continuation in n 


Z'^ — 1 

lim-= log(Z’) = W . 

n^O n 


(1.7) 
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Alternatively, in some sitnations (inclnding all the cases we analyze in this paper), the free 
energy can be snfhciently controlled directly, and then the replica trick is not necessary. 

Analyzing disorder in free held theories (with disorder for the free scalar held) is 
relatively straightforward, and we begin by reviewing what is known about this case. We 
add some new observations on the nature of the low-energy hxed points (when these hxed 
points exist). In particular we show that they are not scale-invariant. 

However, in many experimental situations one is interested in strongly coupled held 
theories, and then there are very few methods for analyzing disorder systematically. For 
example, consider the Ising model in d dimensions with the exchange energy (the bond) 
having a small stochastic component of the type (1.2). This is described in the continuum 
limit by making the coefficient of the energy operator, e(x), a random variable as above. 
For d = 3 this modihes the properties of the phase transition, and one hnds new critical 
exponents. It is of great interest to hnd a useful approach to computing these critical 
exponents. One method that has been attempted is the epsilon expansion around d = 4. 
It turns out that the expansion is in ^/e and it unfortunately does not work as well as for 
the pure Ising model. See [|| for a relatively recent review and references. 

Another method that is useful when the disorder is close to being marginal is an 
expansion in {d/2 — A) (for the Ising model A ~ 1.41, and hence d/2 — A ~ 0.09 is 
numerically rather small and one can hope that expanding in it would give reasonably 
good results) - one can view this as an expansion in the heat-capacity critical exponent. 
When ((i/2 — A) is small the disordered hxed point is sometimes close to the original, pure, 
hxed point. This method was used in two dimensions, see for example and an 

analysis in 2 -|- e dimensions was considered in [H . The application of this approach to the 
random 3d Ising model will appear in a future publication (see m for some preliminary 
details). 

In this paper we study a diherent method for controlling the ehects of disorder in 
interacting theories, which is to take the large-A limit of vector models or of gauge theories. 
In most of this paper we analyze large-A theories where the coupling constant of a ‘single¬ 
trace operator’ is taken to be a random variable. We analyze many aspects directly in held 
theory. In addition, if this theory has a known dual given by a weakly coupled gravitational 
theory on AdSd+i, then W[h,g^,g ‘^,...] of (|1.3| ) is equal to the AdS^+i bulk action with 
specihed boundary conditions, and one can consider its average. If the gravitational theory 
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is weakly curved then it is possible to explicitly study these disordered theories in AdS(i_|_i. 
This is done in section 5 (see |1TT|,[T^ for earlier works, and see also H)l 

In the extreme large-A^ limit, the held theory behaves as a ‘generalized free held 
theory’, and we can exactly solve the path integrals that appeared above, and compute 
the correlation functions of the disordered theory. We hnd again that scale-invariance 
of disordered hxed points is not guaranteed (even though they are hxed points of the 
renormalization group). We argue that when the disorder is marginal, it can sometimes 
be an exactly marginal deformation of the disordered hxed point, and that in some other 
cases the low-energy theory is scale-invariant but has unconventional dependence on the 
disorder. It would be interesting to understand if these non-standard properties of large- 
disordered hxed points can occur also for hnite-A^ theories. We also analyze what happens 
when a double-trace deformation is present in addition to the disorder. 

The extreme large-limit is dual by the AdS/CFT correspondence to a free held 
theory on AdSd+i, and in section 5 we use holography to make detailed computations 
and compare them to the held theory expectations. In particular, we reproduce the beta 
functions that the held theory predicts. 

In sections 5 and 6 we study the leading 1/N corrections to the inhnite N results. 
When the disorder is of order one in the large-A^ limit, these give small corrections that 
become important only at distances that (for marginal disorder) are exponentially large in 
N. However, these ehects are important when the disorder scales as the number of de¬ 
grees of freedom (namely, as N‘^ for gauge theories and as N for vector models). We cannot 
compute the 1/N corrections exactly, but we can analyze them perturbatively in the disor¬ 
der, both holographically (where this involves an interacting held theory on AdS^+i) and 
in held theory, and we do this at leading order. In general the behavior found for inhnite 
N is modihed. At leading order in 1/N we hnd unconventional momentum-dependence 
in correlation functions of the disordered theory, with extra logarithms. Holography may 
allow for an exact analysis of this limit in some cases, but we leave this to future work. 


^ See [r^lh,^| for the holographic analysis of disordered theories of the type that is relevant for 
quantum phase transitions, where disorder breaks the Poincare symmetry. See also, for instance, 


and references therein for additional recent studies of holographic disorder. 
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2. Disorder in Free Field Theories and in Large-Theories 

2.1. A Review of Disorder in Free Field Theories 

We begin by analyzing disorder in free scalar field theories. This can be a first ap¬ 
proximation to disorder in weakly coupled theories. It is also relevant for second order 
phase transitions when we have, on one side of the transition, a broken continuous global 
symmetry, since this gives a free theory of Nambu-Goldstone bosons at low energies. 

Our action is 

S = J d‘^x + h{x)a{x)^ , (2.1) 

and we take the coupling h to vary randomly as in ( |1.2|) . For specific h{x) there is a 
competition between the fact that the field <j{x) wants to align itself with the local source 
/i(x), and the kinetic term that suppresses variations of the field. 

For a given h{x), we can compute the (normalized) two-point function of a{x) in 
closed form,@ 


{a{x)a{y))h = 


1 


{x — y)‘^ ^ 
Thus, we find the exact result 

{a{x)a{y)) = 


+ I d'^zd'^w 


h{z)h{w) 


{x — zY ‘^{y — wY ^ 


(2.3) 


(x - y) 


d-2 


+ / d^Z 


{x — zY “^{y — z) 


d-2 


(2.4) 


The second term on the right-hand side is infrared divergent for d < 4. Therefore, the 
theory (|2.1|) with h{x) a random field is sick for d < 4. This is the simplest example of 
a conformal field theory that ceases to exist because one of its couplings is taken to be 
a random field In particular, this computation shows that phase transitions with 

spontaneous breaking of continuous symmetries cease to exist in d < 4 once we make the 
couplings to the Nambu-Goldstone bosons random 


^ The propagator in position space of a free field normalized as in (2.1) is 

r(d/2-l) 1 


(a(x)o-(y)) = 


47r‘^/2 (a: — yY~‘^ 


( 2 . 2 ) 


We omit these numerical prefactors and some combinatorial coefficients in this section in order 
not to clutter the equations. 
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We can see another manifestation of the sickness encountered in (p.4|) by computing 


= j ■ ( 2 - 5 ) 

Therefore, 

. ( 2 . 6 ) 

This is again infrared divergent for d < 4. 

Note that the connected average {a{x)a{y)) — {a{x)){a{y)) is altogether independent 
of disorder and hnite! This is the combination that arises from the second derivative of 
the disordered free energy with respect to an extra source for a{x). 

However, it is possible to see the infrared catastrophe also at the level of the disordered 
free energy. To that end, we introduce a source for the composite operator cr^(a:), 

S'= y d‘^x ^{da{x))'^ + h{x)a{x) + g{x)a'^{x)^ . (2.7) 


When the CFT describes a second order phase transition, g is proportional to (T — T^). 
Averaging over the random held h and computing the hrst derivative with respect to g of 
the disordered free energy we get 


SWp 

Sg{x) 


((t2(x)) 



h{y)h{y') 
yY~‘^{x — y'Y-"^ 



c 


2 


_ y^2d-4 ■ 


( 2 . 8 ) 


The integral is infrared divergent for d < 4. Hence, the operator a‘^(x) does not exist as 
a local operator in the disordered theory, and the theory is thus sick. This is very similar 
to the familiar infrared divergence that leads to the Mermin-Wagner-Coleman theorem in 


d = 2 [| T ||2 


2.2. The Low-energy Limit for d > 4 

In the cases where the theories described in the previous subsection exist, namely 
d > 4, it is interesting to ask what is their behavior at low energies. 

Consider the connected averaged two-point function 

{aYx)aYy))-{cjYx)){cTYy)) ■ (2.9) 

According to what we found above, this should be well dehned for d > 4. In second- 
order phase transitions we can interpret this correlation function as measuring the second 
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derivative of the disordered free energy with respect to the temperature, namely the heat- 
capacity exponent. One hnds that in the deep infrared (compared to the scale set by c), 
and for d > 4, 


{x)a‘^{y)) - (cT2(x))((T2(y)) ~ 


1 


{x-y) 


d-2 


(TZ 


{x — zY — z) 


d-2 


_ y'j2d—6 


[X 


( 2 . 10 ) 

The scaling dimension of in the infrared is thus shifted (in the second derivative of 
the free energy) from (d — 2) to (d — 3). As a result, the heat-capacity exponent of the 
disordered theory vanishes in d = 6. In d = 6, the two-point function ( |2.1CI|) can be written 
as ~ / d^pe*^^ log(p^). This “dimensionless” behavior in d = 6 is the starting point 

of a systematic e-expansion of the disordered free energy of interacting theories around 
d = 6 ||2^. It is important to point out that again the connected average (|2.10|) behaves 


very differently from {a‘^{x)a‘^{y))^ which is dominated, like the disconnected contribution 
(cr2(a:))(a2(p)), by a term that scales as Y/[x — at large distances. 

Let us also discuss the connected averaged three-point function of namely, 
-TJTe- Qj^g hnds that it receives a contribution that goes like Y from disorder, but 


Sgi3:)Sg(y)Sg(z) 

all the contributions of order Y and higher exactly cancel: 




D 


1 


59{x)5g{y)5g{z) 


[x — y)‘^ “^{x — z)‘^ ‘^{y — zY ^ 
^2 


+ 


[X 


yY ‘^{x — zY ^{y — Y 


d-2 


-|- 2 permutations 


( 2 . 11 ) 


In the infrared the terms proportional to Y always dominate. As expected, for d > 4 
the correlation function ( p.ll|) decays at long distances and so the theory is well behaved. 
However, in the deep infrared the three-point function (|2.11|) is inconsistent with conformal 
invariance. In fact, it is not even consistent with scale-invariance.S Based on the two-point 


^ Note that in momentum space there is a well-known relation between the IR 

limit of disorder-averaged connected correlation functions of the d-dimensional free theory, such as 
( 2.10 ) and (|2.11 ), and connected correlation functions of the pure (d— 2)-dimensional free theory. 
This relation holds when all momenta lie in a (d — 2)-dimensional subspace, and it continues 
to hold to all orders in perturbation theory (though it does not hold non-perturbatively). This 
relation implies that these specific disorder-averaged connected correlation functions in momentum 
space obey a (d — 2)-dimensional scaling symmetry (and even a conformal symmetry), as do the 
correlation functions (|2.10 ) and ( p.ll ) when integrated over two of the dimensions. It also implies 
some relations between the critical exponents of the d-dimensional disordered theory and the 
(d — 2)-dimensional pure theory. However, the full correlation functions of the disordered free 
theory do not obey any scaling symmetry, as exhibited already by the form of (^.1(]|) and (2.11). 











































function ( |2.1C1|) one would have expected that cr^ carries dimension d — 3, but the overall 
dimension of the three-point function (|2.11|) in the infrared is Sd — 8 rather than Sd — 9 i. 

In the next section we will see how these results are reproduced using the replica trick, 
and we will return there to the discussion of scale-invariance in the infrared. Here we only 
mention that this non-scale-invariance in the infrared can be traced in the replica trick 
computation to the fact that the replicated theory is unstable; the deep infrared limit and 
the n —)■ 0 limit of (1^) do not commute. 


2.3. Disorder in the large-N limit 


Large-held theories (such as SU{N) gauge theories in the ’t Hooft limit or 0{N) 
vector models) provide useful examples where disorder can be analyzed explicitly in the 
N ^ oo limit, since connected correlation functions in these theories are suppressed by 
powers of 1/N. 

The simplest case arises when we have disordered couplings for ’single-trace operators’ 
(operators invariant under SU{N) or 0{N) that cannot be written as products of invariant 
operators) in a large-theory and we normalize their two-point functions and sources to 
be of order one in the large-A^ limit; for inhnite N this leads to ‘generalized free helds’. 
Generalized free helds are CFT operators that obey the same statistics as free helds, even 
though they can have diherent dimensions from free helds. We normalize the two-point 
function of a scalar operator 0(x) of dimension A to 


{0{x)0{y)) 


1 

{x-yy^ ’ 


( 2 . 12 ) 


and higher-point functions are then given by Wick’s theorem. The only operators in such 
theories are products of O and of derivatives of O. As a consequence, the OPE of 0{x)0{y) 


® Another option is to force the connected correlation functions to be scale-invariant in the 
infrared, by scaling each operator according to the scaling dimension implied by its two-point 
function. With this scaling we hnd that all connected n-point functions with n > 2 (such as 
( 2.11 )) vanish in the IR. This would imply that the disorder-averaged IR theory is non-local. For 
instance, all connected 3-point functions {a{x)a{y)0{z)) with any operator O would vanish in the 
IR, implying that a{x)a{y) has no overlap with any local operator (except the identity operator); 
this is unlike the situation in the UV theory where the connected {|cr(a:)cr(y)cr^(z)^ is non-zero. 
The same is true also for any other product of two operators 0\{x)02{y) in the IR. Thus, the 
interpretation of having no scale-invariance in the IR seems more reasonable to us. We thank S. 
Rychkov for discussions on this point. 
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does not include the energy-momentum tensor (unless A = (d — 2)/2, in which case O is 
an ordinary free held). Thus, one cannot think of such a theory as a local CFT, but 
one can think of it as a decoupled sector of a local CFT with ct = oo, where ct is the 
central charge appearing in the two-point function of the energy-momentum tensor. In the 
AdS/CFT correspondence, a generalized free held is related to a free held of mass m in an 
AdSd+i space of radius L, with A(A — d) = m^L^.0 

In this section we study such theories with arbitrary values of A < d/2, and with 
disorder for O. In this case we can compute everything exactly as a f un ction of c^, as 
long as remains hnite in the large-A limit. As in the previous subsection, we hnd a 
non-trivial low-energy limit with novel behavior. In the next section we will discuss the 
specihc case of A = d/2, where some new features arise. In large-A disordered theories 
the free energy and correlation functions can be exactly computed as a function of the 
interesting couplings, so we do not need to use the replica trick. However, in the next 
section we will perform computations also using the replica trick, and we will see that in 
some cases there are subtleties in using the replica trick for generalized free helds. 

The computations in this subsection are straightforward generalizations of the compu¬ 
tations performed above for free helds. Consider a generalized free held O with A < d/2, 
so that disorder is a relevant perturbation. Its two-point function in the presence of some 
source h. 


{0{x)0{y))H = 




(2.13) 


can be solved exactly. In momentum space we get (up to unimportant prefactors that we 
will ignore throughout this section) 


{0{k)0{l))h = + /) + k^^-^l^^-‘^h{k)h{l) , (2.14) 


and in position space 


{0{x)0mn 


2A 


-f / d'^zd'^w 


d.d.. h{w)h{z) 


X 


{x — z)‘^^w'^^ 


(2.15) 


We can now directly average these expressions using ( |1.2|) to hnd an exact result [jH 

(^2)d/2-A _|_ ^2 


{0{k)0{-k)) = 




d-2A 


(2.16) 


^ Such relations between sectors of CFTs and held theories on AdS space can exist also for 
non-free held theories on AdSd+i 
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In the deep UV we find 


lim {p{k)0{—k)) = —)■ lim {p{Q)0{x)) — x (2-17) 

k^^oo X^O 


and in the deep infrared 



= c X 


2d-4A 


(2.18) 


We see that in the UV, there is no effect of the weak disorder, and the two-point fnnction 
in the UV is the same as in the pure theory, as expected when the disorder is relevant. In 
the infrared, however, the disorder term dominates and the behavior is altered. 

The distribution is non-integrable for d — 4A > 0, i.e A < d/4. Exactly 

when this happens the correlation function in position space does not decay. This signifies 
a disease of the theory which generalizes the result of ||2^ (our computation reduces to [p0| 
when A — [d — 2)12, as reviewed in subsection 2.1). Note that if the original GET is 
unitary, then A > d/2 — 1, and thus d — 4A > 0 can only happen if d < 4. 

Let us next consider the correlation function (|2.16|) from which we subtract the piece 
which is disconnected before the d-average is taken. Indeed, {0{k))h = k‘^^~'^h{k), so 
{0{k)){0{l)) = Hence, the piece that is dominant in the infrared in (|2.1(j| ) 

exactly cancels with the average of the disconnected correlation function, and the second 
derivative of the free energy behaves as 


{0{k)0{-k)) - {0(k)){0(-k)) = , 


(2.19) 


which is completely independent of disorder! It is easy to see that the same is true also for 
higher derivatives of the free energy, and even for connected correlators of O with other 
operators. This is due to the fact that all the connected correlators of more than two 
generalized free fields vanish. 

As in subsection 2.1, to see the sickness of the theory for A < d/4, we can also add a 
source for the operator O^. Then we can take a derivative of the free energy with respect 
to the source and obtain 



( 2 . 20 ) 


If we average over h, only the zero momentum piece has a non-zero average (as expected 
from translation invariance of the statistical ensemble) and we find 



( 2 . 21 ) 
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This is infrared divergent if A < (i/4. So we conld say that the operator does not exist 
in the disordered theory if A < (i/4. For example, consider the free (non-critical) 0{N) 
vector model at very large N. The operator O = is a, dimension A = d — 2 generalized 
free held operator. Its dimension is smaller than d/4 for d < 8/3. Hence, in this range, 
introducing small disorder in cp would render the theory ill-dehned. 

Now we will assume d/2 > A > d/4 and compute the infrared anomalous dimension 
of O^. We will study the averaged connected two-point function. This subtracts the most 
infrared divergent pieces and has a clear physical interpretation as the second derivative 
of the free energy with respect to some source. The hnal answer for the average of the 
second derivative of the free energy with respect to the source of is 

{0‘^{k)0‘^{-k)) - {0‘^{k)){0^{-k)) - (fc2)-rf/2+2A ^ ^2^^2^-d+3A _ ^2.22) 

There are no higher corrections in c^. The dominant term in the UV is the hrst one which 
leads to the dimension 2A, as appropriate in a generalized free held theory. The infrared 
dominant term is the second term. Thus, behaves as a held with a shifted dimension: 


Ao2 : 


2 A 


3A - d/2 . 


(2.23) 


As expected, the marginal case A = (i/2 is a hxed point of the transformation. Also, the 
dimension of always decreases under such renormalization group (RG) hows. 

However, the infrared theory does not adhere to the rules of any standard scale- 
invariant hxed point. Indeed, if we consider higher derivatives of the disordered free energy 
with respect to the source of they are all at most linear in c^, and we get in the deep 
infrared higher correlation functions of which do not respect the scaling (|2.23| ). For 
example, generalizing (|2.11| ), the averaged connected three-point f un ction scales in the 
infrared as 


{0‘^{x)0‘^{y)0^{z)Y 


{x — '^{y — zY^{x — z) 


\2A 


-|- 2 permutations . (2.24) 


Such a three-point function is certainly inconsistent with being a primary of the confor¬ 
mal group, and, given ( p.23|) , it is even inconsistent with being an eigen-operator of the 
dilatation operator in a scale-invariant held theory. As in the free held theory discussed 
above, one can also interpret the theory as having scale-invariance but with all higher- 
point functions vanishing, but this seems to lead to a loss of locality. We will discuss the 
interpretation of these results in the replica trick computation in the next section. 
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2.4- Generalized Free Fields with A > d/2 

In this case disorder is irrelevant according to the Harris criterion so it would appear 
that there is nothing to be discussed. We will see below that this is not quite the case. 
The momentum space two-point correlation function in the pure theory is given by 

{0{k)0{-k)) = k^^-^ . (2.25) 


However, one can also add a constant which corresponds to a contact term in position 
space 


{Oik)0{-k)) - 


(2.26) 


where Ajjv is some unknown UV scale which is scheme-dependent (for instance, it can 
depend on the details of an underlying lattice theory).! 

If we deform our theory by disorder for O, then, for instance, the disordered two-point 


function (|2.16|) now becomes (up to unimportant numerical coefficients) 


{0{k)0{-k)) = . (2.27) 


Keeping the factors that are most important in the infrared we obtain a contact term 
plus (1 -F -F ■ ■ -. As long as Kuv 7^ 0 the effects of disorder are thus not 

suppressed compared to the pure theory.i On the other hand, if we choose Auv = 0, the 
term proportional to would be multiplied by which is parametrically suppressed 

in the infrared. One can verify that many other correlation functions (including connected 
averaged correlation functions) are affected by disorder, and that the effects of disorder 
cannot be removed by a redehnition of the operators when Ajjv 7^ 0. 

We thus see that although the disorder is irrelevant in the technical sense, it does 
modify correlation functions in the infrared. (This is perhaps reminiscent of dangerously 
irrelevant operators in standard renormalization group flows. The difference is that here 
the effects of irrelevant disorder are only present when we allow for a contact term as 
in (p:^.) 


® If one imagines coupling 0{x) to a source h{x), the above contact term is obtained by adding 
to the Lagrangian ~ A^^“^ J d‘^xh^{x). Such a contact term would be generically induced by an 
RG flow. We will see in section 4 that this indeed happens in the double-trace flow. 

^ In position space, the contribution containing arises from the integrated disorder 

operator coinciding with one of the insertions in the correlation function. 
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3. Disorder at Large-A^ using the Replica Trick 

3.1. Large-N Theories with Disorder for an Operator with A = d/2 

In the marginal case A = d/2 it is natnral to compnte the beta fnnction for the 
disorder. We can compnte this using the replica trick (|1.7|) . The approach we use below 
is analogous to that of The results are however not the same because in the 

context of large-A theories double-trace operators play an important role. 

We start by considering 


W^ = 


r h(x)h(x) 

[Dh]Z’'[h]e~ J * 2 .'‘ , 


(3.1) 


where 


Z[h] = J 

We can then compute the averaged disordered free energy (|1.3|) by 


Wd = 

dn 


n=0 


We write 


/ n 

A=1 


(3.2) 


(3.3) 


(3.4) 


where capital Latin indices go over the replicas, A = 1, • • •, n. We therefore have the path 
integral 


/ n 

A=1 


We solve the h path integral hrst. The equation of motion of h sets h{x) — —c^ YIa ^a{x) 
and thus 

Wn= f (3.6) 

A=1 

Therefore, we have a collection of n identical CFTs perturbed by some marginal couplings 
that connect them. 

In (|3.6|) one has to decide what to do with the operator Oa{x)Oa{x) for some partic¬ 
ular A. In generic held theories there is no marginal operator appearing in this Operator 
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Product Expansion (OPE), so this term can just be dropped. However, in large-theories 
(generalized free held theories) the OPE takes the form 

C>(x)C>(0) ~ ^-fC>2(0) + ... . (3.7) 

Thus, we simply interpret the product Oa{x)Oa{x) as the operator 0\{x). With this 
interpretation we have (keeping only marginal terms) 




f[ [V^iaW ^_ 

A=1 


Expanding (|3.8|) in weak disorder, we get 


(3.8) 


^ A=1 

2 r 4 

l + y / d'^x Oa{x)Ob{x) + y 
d A,B 

(3.9) 

The third term in the second line can renormalize the second term in the second line. In¬ 
deed, marginal terms in the OPE {OaOb) (a^) {OcOd) {y) lead to logarithmic divergences 
at short distances, which can be absorbed either by renormalizing or by adding new 
terms to the action. 

In large-theories there is just one way to get a marginal term as x —)■ ?/ in the third 
term of (|3.9| ), i.e. from the unit operator in a single contraction of (P’s. This gives 

{OaOb){x){OcOd){0) ~ ^ {SacObOd + 3 permutations), (3.10) 


dWy OA{x)OB{x)Oc{y)OD{y) + 


A,B,C,D 


and 


OaOb{x) Y, OcOoiO) ^^aObW . (3.11) 

A,B C,D ^ A,B 


The integral d'^x over (|3.11|) 
have a beta function0 


in (|3.9|) diverges logarithmically, so the conclusion is that we 
dr^ 

=-7(d)cS + 0(c6) , (3.12) 


Similar methods were applied in d = 2 to describe the disordered Ising and Potts models, 
see. e.g. [|,@,|,|,|9l • 
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where 7 (rf) = Vol(^^-^) = 27r'^/Vr(rf/2). 

The beta function for the physical disordered theory is obtained by simply substituting 
n = 0 in ( p.l2|) , letting the derivative in ( |3.3|) act on the hrst line of (|3.9|) , so 


dc2 


(3.13) 


This shows that in generalized free held theories the marginal disorder c is exactly marginal! 

At hrst sight this is not too surprising, since we found in the previous section that 
the averaged connected two-point function of O does not depend on the disorder (so one 
could say that O does not get renormalized as a result of disorder). However, the averaged 
connected correlation f un ction of double-trace operators now takes the form 


{o\x)o^m - mx))mo)) = ^ + f ^7 


y'‘{y - xY 


(3.14) 


Equivalently, in momentum space. 


(C>2(fc)C>2(-fc)) - (C>2(A:))(C>2(-fc)) = 

2 j d‘‘q\og{q^ / iY)\og({k - qf ! + AY j (Y qlog(q^ / Y) \oY‘{{k - qY / Y) . 

(3.15) 

Clearly the hrst term makes sense because if we rescale p we get an integral that does 
not depend on k and hence a contact term. But the second integral transforms under 
rescaling p non-trivially already at separated points. Therefore, the two-point correlation 
function looks reminiscent of a logarithmic CFT (logarithmic CFTs are argued in to 
arise naturally in disordered theories): 


(O2(x)(92(0)) - ((92(x))((92(0)) = 

(92 may seem to have a non-zero anomalous dimension as a fu n ction of the exactly marginal 
disorder c^, but the result (3-16|) is actually exact for any so it does not exponentiate 
as an anomalous dimension should. 

To summarize, the exactly marginal parameter does have an ehect on the correlation 
functions. 
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3.2. Exact Replica Computations and an Interpretation of non-Scale-Invariance 

Next, we want to understand from the point of view of the replica trick why the 
disorder limit n —)■ 0 of (|3.8|) leads to a non-scale-invariant theory for the case of A < d/2. 
We can solve the replicated generalized free held theory (|3.8|) exactly for any n. For any 
integer n 7 ^ 0 we have a standard QFT so we expect the low-energy limit to be scale- 
invariant, but this is not what we found for the disordered theory in the previous section. 
We will see that the limit of n —)■ 0 does not commute with the long distance limit in the 
class of theories we study here .0 

We study the theory of n initially decoupled generalized free held theories with an 
interaction ( — b ^aOb), where A, B = 1, ■ ■ ■ ,n. We will calculate the correlation 
function a ^a{x) b ^B{0)) exactly. By the time we take n —)■ 0, the derivative of 
this with respect to n gives the two-point connected average of the generalized free helds 
in the disordered theory. Note that this correlation function in the replicated theory is not 
normalized, so we always have to multiply the answers below by ( 1 ) at the end of the day. 
That, however, will not change our conclusions. After some combinatorics we hnd 




n 


B 


kd 2A _ 


(3.17) 


This two-point function signals an instability because of the singularity in the denominator. 
This is not surprising as the interaction term that we must add in the replicated theory, 
J2ab ^aOb, appears with the wrong sign in the potential Now: 


A. If we hrst apply to (|3.17|) the operator we clearly get 

1 


{o{k)o{-k)y 


J^d-2A 


(3.18) 


This corresponds to O having scaling dimension A (in agreement with ( 2.19 )). Note 
that it is completely independent of disorder, as we found above. 

B. If, on the other hand, we hrst take the limit of small momentum (crossing the singu¬ 
larity at kd~‘^^ = nc^), we get 


iJ20A{k)J20B(-k)) = 


1 


1 


B 


nc^ 


k 


d-2A 


+ 


(3.19) 


Even for integer n the replicated theory is not really a standard QFT in the case of generalized 
free fields, because of an instability discussed below. 
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The first term is a contact term, and the other terms have singnlar dependence on n 
for small n. This formnla naively implies that the infrared dimension of (9 is d — A, 
which is what we are nsed to from donble-trace deformations (thongh for any n ^ 1 
the deformation is not eqnivalent to a donble-trace deformation). 


We therefore clearly see that the limit of small n and the infrared limit do not com- 
mnte. This is dne to the singnlarity in (|3.17|) , which signals an instability at lower and 
lower energies as n —)■ 0. If we go to very small momenta hrst, we go below the nnphysical 
pole (whose location scales as a positive power of n). Bnt if we take n —)■ 0 hrst, the 
instability is “seqnestered,” and we get the correct answer for the disordered behavior. 

The analysis is particnlarly simple for free scalar helds as in snbsection 2.1. In this 
case the replicated theory is again a free held theory, with one held of mass m? — —nc^, 
and with [n — 1) massless free helds. The presence of the tachyonic held for any hnite n 
indicates that the theory is nnstable, and this becomes problematic once we go to momenta 
below /c^ = nc^. 

We now tnrn to the correlator {J2a After some combinatorics 

one arrives at the expression 


j2ol(k)j2o%(-k)) = 2 j 

A B J \ 


-I- 


+ 2nc\k-l) 

^4A-2d 


2A-d 


4A-2d 


1 — ^ 


(k-l) 


4A-2d 


1 — ^ 1 — nc^{k — 1) 


2A-d 


(3.20) 

Now we demonstrate again that the infrared limit and the small n limit do not commute: 
A. If we apply the operator hrst we get 


{o^{k)o^{-k)) = 2 1 d^^ii p^-^^ik - + 2c\k - ly^-’^i 


2A-d 


\2A-d74A-2d 


(3.21) 


This coincides with (|2.22|) , and the infrared dominant term is the one proportional to 
c^, which gives that translates in position space to This is precisely 

the transformation ( lj.23|) . 

B. If we take the infrared limit hrst, we notice that for instance for n = 1 the leading 


infrared term cancels between the three terms in ( p.20|) . This is important, because in 
this case where we just have an almost standard double-trace deformation, we indeed 
expect that the dimension of should be 2d — 2A. For any other n > 0 the 
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cancellation does not happen and we get that the leading infrared behavior is 
corresponding to x~‘^^ in position space, so that the dimension of stays 2A. 


It wonld be interesting to nnderstand if the breakdown of scale-invariance in disordered 
theories is always associated to snch instabilities in the replicated theory. 

Preliminary stndies snggest that the issues of non-scale-invariance ( |2.24|) and logarith¬ 
mic correlation functions as in ( b.l6|) can also be naturally understood by generalizing the 
Osborn equations to disordered systems. However, we leave this topic for the future. 


4. Double-Trace Deformations and Disorder 


We now generalize the study above to include a double-trace operator 0‘^ in the action 
in addition to the disorder. The reasons for doing that are twofold: 

1. In a general held theory we expect that such double-trace deformations should be 
present in any physical application (since the double-trace deformation is relevant at 
large N if and only if disorder is relevant). 

2. We have seen above that scale non-invariance is related to an instability of the repli¬ 
cated theory. We would like to test this by curing the instability with an explicit 
double-trace deformation. 

Our theory is thus 

S = Sq + -A f d'^xO‘^{x) -\- f df‘xh{x)0{x) , (4-1) 


with h a random Gaussian held. Repeating the steps near ( pj.l6|) we hnd that the two-point 
function now takes the exact form ||11|| 


{0{k)0{-k)) 


(fc2)d/2-A + + A)^ ' 


(4.2) 


For A = 0 it coincides with (|2.16|) , and for c = 0 it coincides with the familiar results in 
theories with double-trace deformations (see e.g. 0). 

We assume A < d/2. At very high energy we always asymptote to (/c^)^“'^/^, which 
corresponds to O having dimension A. This is because the disorder and double-trace 
deformations are both small in the ultraviolet. 
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Let us now compute the averaged connected two-point function. We subtract the 
average of one-point functions from (|4.2|) to hnd 


{0{k)0{-k))-{0ik)){0i-k)) = -^ 


1 


kd-2A ^ 


(4.3) 


So we see that the connected averaged two-point function does not see the disorder at all, 
as in the case of A = 0 ( p.l9|) . 

As before, we can also consider the two-point function of removing the pieces which 
are disconnected before the h-average. One can hnd exact results even in the presence of 
the double-trace coupling. Our hnal answer is 


{o^{k)o^{-k)y 


= / d^i 


1 


+ c^ 


1 


ld-2A ^ X {k - + X \ld-2A^xJ (fc-/)^-2A + A 

(4.4) 

This coincides with the previous result (| 2 . 22 |) for A = 0, where the dimension in the 
infrared is given by 3A — (i/2. Another special case is c = 0. In this case at large 
k, the virtual momentum / is going to be large as well, and we can neglect A and get 


jd^i 


ld — 2A^J^_l'jd—2A ^ —4A 




This power-law behavior is associated to an operator 


of dimension 2A, as expected from at high energies. At low energies all the momenta 


are small and we expand and get J d^l 


1 I 

A 2 A3 


+ 


^d-2A(^_^)d-2A 

A4 


-!-■■■. The hrst two 


terms do not give nontrivial momentum-dependence in the infrared, and so the leading 
piece comes from the third term, which gives —)■ Hence, the infrared 

dimension of the operator is given by 2d — 2 A, as we expect for the ordinary double¬ 
trace deformed theory. 

Analyzing (|4.4D for generic non-zero c and A, it turns out that the infrared scaling is 
always so that has dimension 2d — 2A. This has a simple interpretation. If 

we imagine hrst reaching the hxed point of the double-trace deformed theory, then O has 
dimension d—A. Therefore, according to the Harris criterion, the disorder is now irrelevant, 
and hence the two-point function after turning on disorder has the same IR scaling as in 
the ordinary double-trace deformed theory. Therefore at the level of two-point functions 
the qualitative picture of the renormalization group how is as follows: 
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Disordered fixed point, 
O: A ( 

O^: 3A-d/2 


/ / / ' 
I / / / I 

I I I I 

I y I 

\ ' \ I 

\ \ \ ' 
\ \ \ \ 


\ \ 


UV pure fixed point, 
O: A 
O^: 2A 

0=0 flow 


IR pure fixed point, 

O: d-A 
O^: 2d-2A 

Figure 1: A schematic picture of the large-N renormalization group flows in the 
presence of disorder and a double-trace deformation. 

However, in view of subsection 2.4, this is not quite the whole story. In higher-point 
functions we now hnd results that are consistent with scale-invariance but depend non- 
trivially on the dimensionless parameter c^/A. At long distances we have at separated 
points an exact result 


{Oflx,) • ■ ■ = {Oflxi) • ■ ■ {1 + 2 


X 


(4.5) 


The dependence on c^/A cannot be removed by a redehnition of the operators. From the 
point of view of the infrared double-trace hxed point, A is a cutoff scale in the ultraviolet. 
Even though is technically irrelevant in the double-trace hxed point, the additional 
dependence on c^/A in is present because the RG how to the double-trace hxed point 
induces a contact term as in ( p.27|) . 

Let us note that in the presence of the double-trace deformation, the instability found 
in ( 3.171) is removed. Given any A > 0, for small enough but hnite n, there is no unphysical 
singularity, because the denominator in (p.l7|) is modihed as cn^ —)■ cn^ — A. This is in 
accord with the fact that scale-invariance is not violated by disorder in the presence of A. 


f.l. Marginal Double-Trace Deformations and Disorder 

If we assume that the generalized free held O has A = d/2, then disorder is marginal 
and also the double-trace deformation is marginal. We should thus study the beta function 
in this case, as we did in subsection 3.1. 

In the replica trick, we should now deform the n decoupled theories by 

«S= (4.6) 

A A,B 
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The sign of the double-trace deformation is the usual one necessary for stability if A > 0. 
Simple combinatorics shows that the singular terms in the OPE take the form 

A B A 

^O^^(0)^OcOMx)=^^Zo^Ob(0). (4.7, 

A C,D A,B 

OaOb(O) OcOoix) ~ OaOb(O) . 

A,B C,D A,B 

Combining these, and generalizing the analysis of subsection 3.1, we hud the coupled beta 
functions 


dX 


dlog(/j.) 
1 dc^ 


= 7(d)X^^ 


(4.8) 


2 dlogM 

The physical disordered theory is obtained for n —)■ 0 and we get 

dX 


d\og{n) 

dc^ 


= 7{d)X^ , 

= 2'y{d)c^X . 


(4.9) 


dlog(,u) 

If we take A = 0 initially, then A remains zero and disorder is exactly marginal, as we 
have already seen in subsection 3.1. For any A > 0, A is marginally irrelevant and goes 
like 7 (d) A -A- —1/ log(/r) in the infrared. Plugging this into the running of disorder, we see 
that disorder behaves like 

- —2TT 

log (^) 

in the infrared. Disorder is thus marginally irrelevant in the present situation. 

In order to make contact with the previous subsection, let the operator O now have 
dimension A = in the ultraviolet (e > 0). The beta functions ([4.91) are then 

modihed as 

A\ 

= —eX -|- 7(d)A^ , 

(4.11) 

—ec^ -|- 2'y{d)c^X . 


dX 
dlog((U) 
dc^ 


d\og{ii) 

The coupling A then evolves to its hxed point 7 (d)A* = e, as for = 0, and the beta 


function for disorder becomes 


dc^ 


= ec 


This is equivalent to the physics that we 


dlog(M) 

found in the previous subsection: disorder is irrelevant in the infrared according to the 
Harris bound, since the dimension of O becomes d — A at low energies, and thus 
acquires effectively a negative mass dimension e. 
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5. Holographic Analysis of Disorder 

In the AdS/CFT correspondence |]^ , ^ |33[| , the scalar operator O in the d-dimensional 
CFT maps to a scalar field (/> in a gravitational theory on AdSrf_(_i. This field is weakly 
conpled in the large-A limit of the CFT. In the AdS/CFT correspondence it is natnral to 
normalize the AdS action to have a factor of ct in front of it, so that otherwise this action 
remains fixed in the large-A limit. In this convention the two-point fnnction ((9(9) scales 
as Ct, and thns the sonrce scales as l/yTr and disorder natnrally scales as I/ct- We will 
denote the source appearing in the AdS equations of motion as h and denote its disorder 
averages by noting that compared to our CFT conventions in the previous sections 
c^ ~ c^/ct- On general grounds, because connected correlation functions are suppressed 
by extra powers of ct, the corrections from disorder to correlation functions that do not 
include the operator (9 itself will be suppressed by a power of ct- In particular, this is true 
for the two-point function of the energy-momentum tensor, which will only be corrected 
at order one (compared to the leading term of order ct), and it is also true for partition 
functions on spheres, that were computed in (following |]30||). 

As in previous sections, we discuss Euclidean CFTs.lii We choose a random source 
for an operator 0{x) of dimension A 

h{x) = 0, h{x)h{y) = c^d^^\x — y). (5.1) 

In momentum space we have 

h{ki)h{k2)^C^d^‘^\ki + k2). \ki\<kma.. (5.2) 


where we explicitly put in a UV cutoh kmax to regulate various divergences that we will 
find. Obviously such a cutoh exists in any physical application. We will be interested in 
particular in how the couplings and the disorder depend on this cutoh. Our starting point 
is very similar to that of, for example, [|I3, but with Euclidean disorder (which preserves 
the complete Poincare group ISO{d)) instead of Lorentzian disorder which is independent 
of time.lii 


See footnote 1 for a clarification. 

One can also study holographic disorder using the methods of the holographic renormalization 


group ; this was done for electric disorder in . 
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We begin our analysis with free fields in AdS, corresponding to generalized free fields 
in the CFT (or to infinite ct)- We then add interactions in AdS, so that we include finite- 
ct effects. In principle in holography we should always include (for finite ct) the coupling 
of the field cf) to gravity, but for simplicity we will work in the limit in which the and 
(j)'^ interactions of cj) with itself are much larger than its interactions with gravity or any 
other interactions, and ignore the latter interactions. Gravitational and other interactions 
can be included by similar methods. 

the classical bulk action 


Note that because in the AdS/CFT correspondence |]^ , |^ ,^ 
is the logarithm of the CFT partition function (at leading order in 1/A^), in the holographic 
description we can compute the connected correlators for quenched random disorder (such 
as (|1.6|) ) simply by putting random sources into the classical equations of motion and av¬ 
eraging over them with a Gaussian weight. Therefore, the AdS/CFT correspondence maps 
the problem of disorder in strongly coupled large-N conformal field theories to standard 
stochastic differential equations. 

Note that at linear order in the equations of motion of any field, it does not matter 
whether we first solve its equations or first average over the disorder; at higher orders this 
does matter, and we should first solve the equations of motion for general sources and only 
then average over the disorder. 

We will take the metric on AdS space to be 


= —^{dz^ -I- dmndx'^dx'^) 


(5.3) 


The boundary of space is at z = 0, while large z describes the IR region. In our approxi¬ 
mation where we neglect the coupling to gravity, this will remain the exact metric. For a 
free scalar of mass m on AdS the dimension of the dual operator is A = | — z^, where 


w\ = 


d2 




(5.4) 


When u < 0(A > ^) the random source is an irrelevant deformation. In order to introduce 
a relevant deformation we have to use the alternative quantization for scalars in AdS space 
so that n > 0 ||^, consistent with unitarity for operators of dimension | — 1 < A < |. 
Note that the unitarity condition imposes u < 1. 

On dimensional grounds we expect that if there is any part of the solution that is 
independent of the cutoffs then it should depend only on cz^ (after averaging over the 
disorder). The perturbative expansion in F will thus break down (or acquire a cutoff 
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dependence) for large ^ when > 0, and for small ^ when < 0, as expected from the 
(ir)relevance of the disorder. 

The linear solntion which is regnlar in the IR for an operator of dimension A with a 
sonrce h{k) in momentnm space is 


= j , (5.5) 

where is a modihed Bessel function, up to a constant for u ^ Q, and (—1) for 

z^ = 0, which we will ignore. 

Near the boundary, for z -C 1/kmax: we expect the solution to be governed by the 
source at the corresponding point, so (j){x,z) should go as h{x)z‘^~^ = h{x)z'^/‘^^'^, and 
this indeed follows from the behavior of the Bessel function at small arguments. More 
precisely, this is true for zv < 0, while for z^ > 0 the solution is dominated by the vacuum 
expectation value (VEV) of the corresponding operator so it scales as For the 

special case of z^ = 0 the leading solution ( b.5|) in this regime goes as h{x)z‘^^‘^ \og{zkmax)^ 
again as expected from the form of the non-normalizable solution for the scalar held. As 
we go into the bulk the source gets smeared. 

In the absence of interactions, the scalar held cj) vanishes after averaging over the 
disorder, so let us compute the disorder average of cp: 


(f)‘^{x, z) 



d^k 

(27r)'^ 


\k\-^’^z‘^K^{\k\zf. 


(5.6) 


If we change variables k^ — u^/z, then the average becomes 


(j)^{x, z) — z^’^c^ 


d'^u 

(27r)^ 


u\-‘^''K,{\u\Y. 


(5.7) 


The integral is hnite for d > 4zz, or A > <7/4. This coincides with the held theory 
result described after (|2.18|) , though we are not computing the same observable here. For 
d = 4z/ there are logarithmic divergences from the |rz| —)■ 0 limit. For v < 1/2 we have 
hnite results for d>2. For z <C 1/kmax the average of goes as z^’^{zkmax)^~‘^’^ when 
z/ > 0, as z‘^'''{zkmaxY z/ < 0, and as c^{zkmaxY^o^{zkmax) for z^ = 0. In this 
region the scalings are just the squares of the scaling of cj). In the region z ^ 1/kmax which 
describes the disordered held theory, the average goes as in all cases, including zz = 0, 
independently of the cutoh. 
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If the field (j) is free in the bulk, then this is the whole story, and in the marginal case 
1 ^ = 0 the disorder deformation seems to be exactly marginal. Since there is no dependence 
on the cutoff, it seems that the held theory stays disordered with the same coefficient at 
all relevant scales. This is expected since the large-beta function for the disorder (|3.13| ) 
vanishes in this case. However, does not have any direct held theory interpretation, so 
in the rest of this section we will compute physical observables to understand better the 
dependence on c^. 


5.1. Generalized Free Fields and Double-Trace Operators 


Let us start from the case of a free held theory in the bulk, in which (|5.5|) is the exact 
expression for the held in the bulk. This case was also analyzed using the replica trick 
n|, with similar results. Up to constants, the expansion of (|5.5|) in moment um space 


m 


takes the form (for A ^ d/2) 


z) ~ z^ ^h{k) -t- z^j3{k) -|- higher orders in (5.8) 

where /3(/c) = \k\'^^~'^h{k) (up to constants mentioned above). /3(/c) may be identihed with 
the VEV of the operator, so that we have {0{k)) = \k\'^^~‘^h{k) for any source h{k) (this 
is determined by dimensional analysis; note that the vacuum expectation value here is a 
held theory vacuum expectation value, not a disorder average). By taking a derivative 
with respect to the source we obtain 


{0{k)0{-k)) = (5.9) 

for any source h{k). This is independent of the source, so obviously we obtain the same 
answer also after the disorder averaging. General n-point functions of 0(x) in this case 
are given by products of 2-point functions, so they are all independent of the disorder; this 
trivially reproduces the held theory result (|2.19| ) (note that the correlation functions we 
compute in this section are always connected). 

In the special case of A = (i/2 we have instead 

z) — z'^^‘^{h{x) log(z/u) -1- /3(a;)), (5.10) 


with a dependence of the VEV on the arbitrary scale /i which must be introduced in this 
case. The form of the solution above implies that, up to a rescaling of 

{o{k))^m = ^og{k/pi)h{k) ( 5 . 11 ) 

By an abuse of notation, from now on we will use k both for the momentum vector and for 
its absolute value. We hope this will not cause any confusions. 
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and 


{0(k)0{-k)) = log(fc/f<) , 


(5.12) 


independently of the disorder, as expected. 

Next, let us add to the held theory action a double-trace interaction S = ^ J d‘^xO^{x); 
if we make A space-dependent then this will also be useful for computing correlation func¬ 
tions of 0^{x), but for now let us just take A to be a constant. Again, A is related to A 
in our CFT analysis by a factor of ct- The analysis will again be very similar to the held 
theory, so let us focus on the most interesting case of A = (i/2 where we expect to see 
some logarithmic running Since A and run, we need to hx some UV cutoh kmax 

as before, where we turn on the disorder with a value c^ikmax) = Cg, and also set an initial 
boundary condition for the RG how X{kmax) = ^o- 

In the presence of the double-trace deformation, (^.10|) is modihed to 


(j){x, z) — z^^‘^{[Xl3{x) + h{x)] \og{ziJ,) -t- I3{x)) -|- higher orders in 


(5.13) 


Using (|5.11|) for the new source in (|5.13|), we have 


P{k) = log 



(X^{k) + h{k)) 




m = {om 


h{k) \og{k/n) 
1 -Xlog{k/n) 


(5.14) 


This determines the 2-point function and all correlation functions of 0{x), in agreement 
with the held theory analysis of section 4. In particular, since the 2-point function is the 
same as the held theory result, and is independent of the source, we hnd the same running 
of A as we did in the held theory, namely that A(^) oc —l/log(^). Note that in general, A 
could also start depending on the momentum due to higher derivative double-trace terms 
that could be generated, but we ignore this. 

If we now compute ( ^.13|) at two diherent scales, one of them being the scale kmax 
where we impose our initial boundary conditions, we obtain (up to higher orders in z) 

(f){x,z) ^ z‘^/‘^{[X{n)/3{fi,x)+ h{n,x)]log{zfi) +/3{fi,x)) ^ 

^ (5.15) 

— Z^ {\XQf3{kmax-i hikmaxi Xy\\og(^zkmax^ Pikmaxi . 


Assuming that we take momenta much smaller than we can use (|5.14|) to hnd 


h{^,k) 


h{kr 


1 - A(^) log(/c/^) 
1 Aq ^Og{k/kmax') 


(5.16) 
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Squaring this, and averaging over the disorder letting it depend on the scale and on the 
momentum (unlike the UV disorder Cq), we obtain 


c^(/r, k) = cf 


1 - A(/r) logjk/^) 
1 -^0 ^^^{k/kmax) 


(5.17) 


We are interested in the /c-independent part of this, corresponding to the Gaussian part of 
the disorder, and indeed for /c -C /r < kmax: the result is independent of /c, and we hnd 






X2 


(5.18) 


When joined with the running of A that we found above, this implies c'^(^) — l/log^(/u), 
precisely as in the held theory analysis ( [1.10|) . Note that for A = 0 we hnd no running of 
c^, as expected. 

In order to compute correlation functions of 0^{x) as in the previous section we need 
to take A to depend on x and to take derivatives with respect to it. The general principles 
of the AdS/CFT correspondence imply that this must give the same answers as the held 
theory analysis above, but let us show this explicitly for the one-point function of 0^(0) 
as a function of the source h{x). To obtain this we hrst compute which 

we can obtain by putting A(a;) = Ai(j(a;) (for simplicity we are not including a constant 
double-trace coupling here). With this coupling, ( |5.13|) becomes in momentum space (at 
leading order in z) 


4>{k, z) — z 


- ydl2 


Ai / + h{k) 


log{zfi) -A /3{k) 


(5.19) 


Using the generalization of (|5.11|) we obtain 


P{k) = log{k/fx) 


Ai / d^ip{l)+h{k) 


(5.20) 


Integrating this we obtain 


d'‘kl3{k) =Ai / <fk\og(klti) / d‘‘W{l) + / d'‘k\og{k/ti}h{k), 


( 5 . 21 ) 


leading to 


d^k(3{k) 


f d‘^klog(k//ii)h(k) 

1 — Ai / d^k log(/c//i) 


(5.22) 
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Plugging this back into (5.2C), we find 


{o(k)) = m = iog(i/M) 


1 — Ai / dH\og{l/ii) 


And, taking the derivative with respect to Ai and setting Ai = 0, we obtain 


{0{k)0^{x = 0)) = log(A;/^) / d^/log(//^)h(/). 


(5.23) 


(5.24) 


This is also the derivative of (0‘^{x = 0)) with respect to h{k) (assuming we only keep 
connected pieces), so we have 


(C>^(a: = 0)) = d‘^k\og{k/iJ,)h{k)j . (5.25) 

This is precisely the expected answer that we would get by contracting the two (P’s in this 
operator with additional O’s, so the results manifestly reproduce the field theory results 
for generalized free fields. It is clear that this will be the case also for higher correlation 
functions, and for theories with non-zero A. 


5.2. Adding Bulk Interactions - General 

Next, let us add interactions in the bulk. If we have some bulk potential (in addition 
to the mass term) V{4>) = Xs4>^/3 -|- A4(/)^/4, then in momentum space (projecting the 
equation of motion of (/> on a specific value of the momentum k) we have 

z^dlcf) + (1 — d)zdz4> — z^k‘^(f) — mALAcf) — Aa^^ -|- . (5.26) 

We will mostly specialize to the marginal case where !S. = dj^ and mAlA = —dA/A. 

The general solution to the homogeneous equation in the marginal case is 

^{k, z) = z^l‘^{AK^{kz) + Bloikz)) (5.27) 


for some constants A and B, while the general solution to the inhomogeneous equation 
with a function f{z) on the right-hand side of (|5.26D (for a specific momentum k) is given 
by 


(j){k,z) 




K, 


\{kz) [ - '^-w f{w)dw — {Iq Kq) 

.L w 


(5.28) 
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plus a homogeneous solution (|5.27|) . Assuming that both integrals converge when ^ ^ oo, 
and that the second integral vanishes in this limit, the IR boundary condition sets 5 = 0. 
And, assuming that both integrals converge when ^ ^ 0 (which will be the case for the 
functions / that we will encounter), the fact that we do not want to have an extra source 
in the UV (which we could swallow into a redehnition of the original source) means that 
to get the precise solution we should choose A such that the solution becomes 


(^(k, z) — z 


- yd/2 


-Koikz) f -^^^w-‘^/^f{w)dw-Io{kz) [ -^!^w-‘^/^f{w)dw 

Jo w w 

(5.29) 

The new contribution to the vacuum expectation value of O is then given by 


{0{k)) = - [ 

Jn W 


(5.30) 


For general values of the dimension of O, (|5.29|) takes exactly the same form but with 
the index of the Bessel functions modihed from 0 to (—z^), and (|5.30|) is replaced by 


{0{k)} = -k-" f 


(5.31) 


5.3. (f)'^ Interactions 

Let us start from the case A 3 
and is given by 


0. In this case the source in (5.26) starts at order hJ, 


f{k,z) = \ 4 ^(jA{k,z) = A 4 


■ d^ki 
(27r)^/^ 

r d'^k 


h{ki)z^'^/‘^Ko{kiz)- 


(27r)'^/^ 


h{k 2 )h{k -ki- k 2 )KQ{k 2 z)KQ{\ki + k 2 - k\z). 

(5.32) 

The solution at order will be given by (|5.29|) . 

If we are just interested in computing {0(k)) at order h^, and then using this to 
compute the disorder average of the two-point function {0{k)0{—k)) at order c^, we can 
average already at this stage the disorder over two of the sources in (|5.32|) . Thus, up to 
this order we can replace (5.32) by 


/(/c, z) = A 4 (/)^(/c, z) = C^A 4 


d^kiz^^/^Ko{kizfh{k)KQ{kz) 


I'zkr, 


(5.33) 


X4h{k)KQ{kz) 


d'^uiz'^/^KoiniY^ 


30 

























(up to a multiplicative constant). 

Plugging this into (|5.29|) , we find that the partly-averaged (j) at order which we 
will denote by (^>( 3 ), is given by 


(/)(3)(/c, 2 ;) =c^X4h{k)z^^‘^ 


-Ko{kz) 

-Io{kz) 


Io{kw) 


owkr, 


W 


Ko{kw) 


w 


d‘^uiKo{ui)‘^KQ{kw)dw- 


d‘^uiKo{ui)'^ Ko{kw)dw 


(5.34) 

Note that the integrals over w indeed converge for finite cntoff, as is necessary for the 
validity of this specific solntion. 


Alternatively, we can compnte the average vacuum expectation valne {0{k)0{—k)) 
directly by taking a derivative of (|5.30|) , and we find 


{0{k)0{-k)) = -c^A 4 


Kpjkwf 

w 


(*wkr\ 


d‘^uiKo{ui)'^dw. 


(5.35) 


Since the integral over ui converges as kmax —)■ C) 0 , it may seem that there is no cutoff 
dependence here (and thns also no interesting momentnm dependence, given that the left- 
hand side is dimensionless). However, this is too naive. The integral over w has two 
regions. When w 3> Ijkmax, the integral over ui is approximately given by the integral 
from 0 to 00 , which is some constant Cq independent of w. However, when w -C 1/kmax, 
the integrand can be approximated by fog^(ui), snch that the integral over ui behaves as 
{wkmax)‘^^og'^{wkmax)- Thus, we have approximately, for k ^ kmax (up to constants) 


{0{k)0{-k)) ~ -CW 4 


w 


(wkmax) log {wkmax)dw 


+ 




'l/fcn 


W 


(5.36) 


C. 


~ \og^ {k/kmax) + O {c^ X^log^ {k / kmax)) ■ 

Note that this result is dominated by the contribntion from w S> 1/kmax which is above 
the effective short-distance cntoff on the radial direction. We will try to understand this 
snrprising resnlt from the field theory point of view in the next section. 

Repeating exactly the same procedure for general values of the operator dimension 
A < (i/ 2 , (|5.35|) is replaced by 


{0{k)0{-k)) = -c^X^k 


-2u 


K^{kw)‘^ 2 


w 


K^{ui)^dw. (5.37) 


w 
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Using the expansion of the Bessel function near a; = 0, K^{x) oc x~'^, we hnd that the 
contribution from the region w > 1/kmax is now proportional to \og{k/ kmax) , 

while the region w < Ijkmax has a contribution proportional to k~^'^. Thus, in this case 
we have 

{0{k)0{-k)) ~ -c^KC^k-'^'^ \og{k/krr^ax) (5.38) 

for some constant , with an extra logarithm on top of the power expected on dimensional 
grounds. Again, we will understand the origin of this logarithm in the next section. 


5.A Interactions 

Now let us take a non-zero value of A 3 , so that the operator O has a non-zero 3-point 
function proportional to A 3 (for simplicity we take A 4 = 0 here). In this case we have a 
source for cj) already at order given by 

f{z) = \2,4P‘{k,z) = X^ j ^^l^^ h{ki)z'^KQ{kiz)h{k - ki)Ko{\k - ki\z). (5.39) 


Performing a disorder average here gives something proportional to 5^‘^^(/c), related to a 
constant VEV for the disorder-averaged (fP' (El- 

Anyway, to go to higher orders we now cannot directly average over this, and we have 
to plug this full expression into (|5.29|) , leading to (up to a multiplicative constant) 


z) = j d‘^kih{ki)h{k - ki)- 

-Ko{kz) [ w‘^^‘^Ko{kiw)Ko{\k - ki\w)dw- (5.40) 

I Jo 

Io{kz) [ — w‘^^‘^Ko{kiw)Ko{\k — ki\w)dw . 

J. w 


At order the source is now given by this solution times the hrst order solution, 
namely 


f{z) = Xz{ 4 P{k,z)){^) = 2 Xlz'^ J d'^k2h{k - k2)Ko{\k - k2\z) J d‘^kih{ki)h{k2 - ki)- 

—Ko{k2z) [ Ko{kiw)Ko{\k2 — ki\w)dw— 

Jo w 

Io{k2z) [ — Ko{kiw)Ko(\k2 — ki\w)dw . 

Jz W 

(5.41) 
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Once again, for the purpose of computing the disorder average of {0{k)0{—k)) to 
order 7?, we are allowed to average over two of the disorder coefficients appearing here. 
There is one average here that gives a divergence, which is when we average over the last 
two sources such that /c 2 = 0. The origin of this divergence is confusing at hrst sight, since 
it does not go away even when we put both UV and IR cutoffs on the disorder distribution 
( 0 ). However, if we follow its origin from the bulk point of view, we see that it arises 
from an exchange of a bulk scalar with spatial momentum /c 2 = 0. Thus, this is an IR 
divergence that will disappear in the presence of any IR cutoff on the full theory. We will 
therefore subtract it in our computation. 

Considering the other contributions to the average, we obtain 


f{z) = Xsicp'^ik, z))(^s) = ‘2c^Xlz'^h{k) I d'^kiKo{kiz) 


—Ko{\k + ki\z) [ o{\k + ki\w) 

Jo w 

Io{\k + ki\z) [ + 


w 

(5.42) 


Next we need to plug this into (|5.28|) to obtain the third order solution, but for us it 
is enough to obtain the VEV of O to third order, using (|5.3U| ). Using this to compute the 
2 -point function, we obtain 


{0{k)0{-k)) = -2 c2 A^ f f d^kiKoikiv)- 

Jo Jo 

—i^o(|^ + ^ik) [ — Ko{kiw)Ko{kw)dw— 

Jo w 

Io{\k + ki\v) [ —— ^-w^^‘^KQ{kiw)Ko{kw)dw . 

Jv W 

(5.43) 

The integrals here converge separately, but as in the previous subsection we need to analyze 
them more carefully, and such a careful analysis again gives a \og^ (k/kmax) behavior. 

Note that in this case, and also in the previous case at higher orders in it is 
important that we hrst solve the equations for arbitrary sources, and only then average 
over the disorder. We would get different answers if we would hrst average and then plug 
the results into the higher order computations. Thus, it is not clear if there is any meaning 
to an ehective disordered solution for the bulk held ^(x, z). 
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6. 1/N Corrections 

In subsections 5.3 and 5.4 we saw how we could perform holographic computations for 
disorder with of order the central charge ct, at leading order in the disorder = c^/ct 
( in principle it is straightforward to extend this computation to arbitrary orders). We 
found several interesting results for this scaling of c^, and in this section we would like to 
understand these results from the held theory point of view. So, we consider the leading 
correction as a fu n ction of c^/ct directly in the held theory. For hnite in the large- 
limit this gives the hrst 1/N correction, that becomes important (in the marginal case) at 
distances exponentially large in N, while for of order ct this correction is signihcant at 
distances of order one, and can lead to changes even if (?/ct is small compared to other 
scales. Below we will present two derivations of our results in parallel. One approach is 
direct, and the other is through the replica trick. It seems like the latter approach is easier 
to generalize to higher orders, while the direct approach removes doubts about the validity 
of the replica trick in this situation. 

In principle it is possible to compute {OO) in a perturbative expansion in the disorder, 
using conformal perturbation theory. At leading order this formally gives 

(0(x)0(!,)> = // + ic" J rfw {0{x)0(y)0(w)0{w)) + ■■■. (6.1) 

The limit when the two points near w come together is singular, but we can introduce 
some UV cutoff on the distance between the two points and use the OPE expansion in the 
original theory near w, which has the general schematic form 

0{w)0{w -|- e) ~ coopP{'w)e^^~'^^. (6.2) 

p 

The identity operator can be dropped from ( |0|) since it gives a disconnected contribution 
already for hxed sources. We will normalize the operators P to have a 2-point function 
with coefficient one. Terms with Ap < 2A are UV-singular, and will depend on a UV 
cutoff n on the distance 1/e. Terms with Ap ~ 2A lead to universal contributions and 
should be kept, while terms with Ap > 2A will be negligible as the cutoff is removed. As 
we are studying large-V theories, we will assume that an operator with Ap ~ 2A exists. 
We denote P = for this particular operator (keeping in mind that at hnite but large 
N the dimension is no longer exactly additive, A 02 — 2Ao ~ c^^). There could be other 
operators of dimension 2A, in particular, single-trace operators. These would modify the 
results below, but it is clear how to generalize our discussion when necessary. 
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The contribution from at the order is 


{0{x)0{y)) = -c^cooo^ / d'^w {0{x)0{y)0'^{w)) = 


1 2 2 2A —A„2 / jd 

= 2 ^ 02 dw- 


(6.3) 


(x — ^o2 (^x — w)^o2 (^y —'w)^o2 


The above result arises in the replica trick as follows. We study n copies of the CFT, 
perturbed by b ^aOb- The terms with A = B are interpreted as the operators 

^ 2 A-A(p^^ To compute the average ( |0|) we then consider the correlator (OiOi) and in 
the end of the computation we take n —)■ 0 (without taking a derivative with respect to n). 
Clearly the only contribution at order arises from J d'^w{Oi{x)Oi{w)Oi{y)), and it is 
independent of n. Therefore, one hnds precisely (|6.3|) . 

The integral over w ( |6.3|) diverges at large w if A 02 < dj^. In the large-A limit this 
happens when A < (i/4, giving another exposition of the problems in that case (that we 
have already discussed at length in previous sections). 

As long as I < A 02 < d the integral converges both for large w and for tc ^ x and 
w ^ y. We get (up to constants) the following leading order in result 


{0{x)0{y)) 


1 (-2 
^ 2 2A-A^2 ^0002 

2 [x — y')‘^^^^o2-<^ 


(6.4) 


This result is to be trusted only for A 02 equal to 2A up to inhnitesimally small corrections. 
It is easy to see both in the direct approach and in the replica trick, that in the absence 
of such an operator there is no contribution at order (? that is independent of the cutoff. 

If the dimension IS.Q 1 is exactly 2A, (H) just gives a hnite contribution going as 
(?(?QQQi l{x — In the large-A limit, when we normalize the 2-point f un ction of 

to one, cooo^ behaves at large-A as \/2(l -|- A/ct + ■■■) for some constant A independent 
of Ct- For ~ 1 only the hrst term remains, while for ~ Ct the second term also 

remains. If << ct but (? jcx stays constant in the large N limit, then we can trust 

our approximation and neglect higher order corrections in c and higher order corrections 


m c 


-1 


We are interested in comparing with the holographic computation which sees only 
the disorder-averaged free energy, so we should subtract from (|6.4|) the disconnected piece, 
which is always present and goes (at order c^) as 


{0{x)){0{y)) 


d’^w 


(x — wY^{y — w) 


2A • 


(6.5) 
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At leading order in l/cr this exactly cancels the leading term above, as in our general 
discussion of generalized free held theories in section 4.0 All remaining terms are of order 
(? jcTi ensuring a hnite limit when (? oc ct with (? jcT 1 . 

Thus, for A 02 = 2 A < d, the connected two-point function at order goes as 

A(? 

CTix-yY^-'^' 

We will now include the effects that are associated to A 02 deviating from 2A. In 
general IS.Q1 = 2 A — B/ct + ■ ■ ■ for some constant B, which stays hnite as ct — t 00. 
In such a case we can write = 1 + (B/ct) log(/i) -!-■■■. The term of order one 

contributes to our computation a term proportional to A as in (| 6 . 6 |) . The S-term gives a 
contribution to the 2-point function (|6.4|) proportional to 

c'^Blog{n\x-y\) 

Ct{x — y)4A-d 

Note that in the holographic computation B should have contributions proportional both 
to A 4 and to A 3 , and in particular we found a contribution of this form proportional to A 4 


in our computation with general dimensions A (|5.38| ). 

Our computation up to now was done for the relevant case, A < d/2. In the marginal 
case A = d/2 we need to reconsider the integral over w that we did above, since it has a 
logarithmic UV divergence as re —)■ x and w ^ y. Note that this divergence has a diherent 
origin from the one leading to (|6.7|) - it is related to the fact that A 02 is close to d (but it 
is still also close to 2A). Note also that the disconnected contribution (|6.5| ) has a similar 
logarithmic divergence whenever A is close to d/2. We therefore need to combine all these 
logarithms carefully. 

In a general large-A theory, the dimension of the single-trace operator O in the large- 
N limit goes as A = d/2 — D/ct + ■■ ■ for some constant D. If we now expand ( |t).5|) in 
1/ct we get (keeping only the leading logarithm) 

_ r . ^2 

(l-f (d - 2A)(log(//(x - w)) + log{^{y - w))) + • • •) 


{0{x)) {Oiy)) = / d^w 


{x — wY{y — wY 


= 3(d-2A) 


l{dy 

(x - yY 


\og\^{x - y)) + O ( 4- 


( 6 . 8 ) 


This cancellation is very easy to see in terms of the replica trick. The connected averaged 
correlation fnnction corresponds to the derivative at n = 0 of '^b which goes like 

in generalized free field theory. This is why the connected averaged correlation functions of 
the single-trace operators do not depend on disorder in a generalized free field theory. 
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dropping the leading term (which eventnally cancels) in the second line. Similarly, the 
integral over w appearing in the second line of (|6.3|) for A 02 ~ d gives (again keeping only 
the leading logarithm) 


3((i — Aq 2 ) 


7 (^) 

{x - vY 


log^(/i(x -y)) + 0 


(6.9) 


so that the total contribntion to {0{x)0{y)) goes as 


(3d-3Ao2 -2(2A-A 02 )) 


iog^(^(^ -y))+o{^ 


( 6 . 10 ) 


{x - vY 

When we snbtract (| 6 . 8 |) from (|6.10|) to obtain the connected two-point fnnction, the con¬ 
tribntion is of order c^/cr and the leading log term takes the form 


(2 A — A 02 ) 


7(d)c^ log^{n{x-y)) 

{x-yY 


( 6 . 11 ) 


This is exactly the same form that we found in the previous section, since the Fourier 
transform of \og'^{n{x — y))/{x — yY is proportional to log^(/c//r). So we identify the 
coefficient we found there as the correction at order 1/ct to the dimension of compared 
to twice the dimension of O. It is indeed well-known from the early days of the AdS/CFT 
correspondence 0 that the diagrams that appeared in the previous section (at order A 4 
and at order A 3 ) contribute to this ‘anomalous dimension’ .0 As in previous cases, the 
two-point function of the form ( | 6 . 11 |) does not look like an anomalous dimension for O 
in the disordered theory. So the interpretation of the computations above in terms of 
anomalous dimensions is not straightforward. To learn about the low-energy behavior we 
need to go to higher orders, and to resum the perturbative expansion in c^/cr. This is 
beyond the scope of this paper, but it would be very nice to do so in the future. 
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